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It has been argued that the existence of old neutron stars excludes the possibility of non- 
annihilating light bosonic dark matter, such as that arising in asymmetric dark matter scenarios. 
If non-annihilating dark matter is captured by neutron stars, the density will eventually become 
sufficient for black hole formation. However, the dynamics of collapse is highly sensitive to dark- 
matter self-interactions. Repulsive self-interactions, even if extremely weak, can prevent black hole 
formation. We argue that self-interactions will necessarily be present, and estimate their strength in 
representative models. We also consider co-annihilation of dark matter with nucleons, which arises 
naturally in many asymmetric dark matter models, and which again acts to prevent black hole 
formation. We demonstrate how the excluded region of the dark-matter parameter space shrinks as 
the strength of such interactions is increased, and conclude that neutron star observations do not 
exclude most realistic bosonic asymmetric dark matter models. 



I. INTRODUCTION 

The similarity of the ordinary matter and dark mat- 
ter (DM) relic abundances suggests a common origin 
for both. This motivates the asymmetric dark matter 
(ADM) scenario, which links dynamically the two mat- 
ter components of the universe. In this scenario the relic 
DM abundance is due to a particle-antiparticle asymme- 
try, in analogy with the relic ordinary matter abundance. 
The DM asymmetry is maintained today due to a parti- 
cle number symmetry governing the low-energy physics 
of the dark sector, a dark baryon number. The excess 
of dark particles over antiparticles may be related to the 
excess of ordinary particles over antiparticles via various 
mechanisms, thus establishing a tight relation between 
the dark and ordinary relic densities (see e.g. [1-3], and 
for a review see Ref. [4]). Moreover, ADM models typi- 
cally predict DM masses of a few GeV. This mass range 
is currently favoured by the DM direct-detection experi- 
ments DAMA, CoGeNT and CRESST [5]. 

Recently, observations of stellar objects have been em- 
ployed to constrain the properties of ADM [6-11]. If 
DM interacts with nuclear matter, then it can scatter 
on the nucleons in compact objects, lose energy, and get 
captured in their interior [12, 13]. The DM capture in 
compact objects may have observable consequences. For 
example, if DM self-annihilates, its capture and annihi- 
lation inside the Sun and the Earth can produce poten- 
tially observable neutrino signals (see e.g. [14]), while its 
capture and annihilation in neutron stars can alter the 
thermal evolution of the latter [13, 15]. However, in the 
ADM scenario, today's universe contains no dark anti- 
baryons, all of which annihilated with dark baryons at 
early times. As a result, no DM self-annihilations can 
take place today. 



Non-annihilating DM captured in compact objects ac- 
cumulates over time, without its density being capped 
by self-annihilations. Depending on the local DM den- 
sity, the accretion rate of DM in a star may be large 
enough such that the DM in the core of the star con- 
denses collectively in its ground state. Fermionic matter 
in its ground state is supported against gravitational col- 
lapse by the Pauli exclusion principle, if the total num- 
ber of particles does not exceed the Chandrashekhar limit 

-^Cha ~ (-^pi/ to ) 3 j where m is the particle mass. Bosonic 
matter in its ground state is supported by the uncer- 
tainty principle, which is overcome if the number of par- 
ticles exceeds the much lower limit A£ ha w (Mpi/m) 2 . 
However, this limit is dramatically increased if repulsive 
self-interactions are present. When the Chandrashekhar 
limit is exceeded, gravitational collapse ensues, and a 
black hole can potentially form, unless some other sta- 
bilising mechanism takes over. This may be the case if, 
for example, the particles under consideration are com- 
posite rather than fundamental. Then the degeneracy 
pressure of the constituent degrees of freedoms can even- 
tually withhold gravitational collapse, even when a much 
higher density is reached. 

Following such considerations, Kouvaris and Tinyakov 
showed that neutron star observations exclude funda- 
mental bosonic ADM having weak-scale interactions with 
ordinary matter but vanishing self-interactions, in the 
mass range 2 keV - 16 GeV [8] (see also [10]). In this mass 
range, a black hole would be expected to form in the core 
of neutron stars located in the Milky Way within their 
observed lifetime. The black hole would subsequently 
consume the stars. Importantly, the excluded DM mass 
range includes the range favoured by DM direct-detection 
experiments. In the case of self-interacting fundamental 
bosonic ADM, Ref. [8] derived constraints on the DM 
self-interaction. 1 
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1 For potential implications of ADM capture in the Sun, see [16]. 
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FIG. 1. DM-nucleon scattering and DM self-scattering in- 
duced by the DM-nucleon effective operator. 

In this paper, we point out that: 

(i) The existence of DM-nucleon interaction, necessary 
for the capture of DM in stellar objects, implies 
the existence of DM self-interactions which, if re- 
pulsive, are in most cases sufficiently strong to un- 
exclude the DM mass range of interest. (If the self- 
interaction is attractive, a different set of consider- 
ations is due.) 

(ii) In many ADM models, DM can co-annihilate with 
nucleons or leptons via the couplings introduced to 
relate the matter-antimatter asymmetries of the or- 
dinary and the dark sectors in the early universe. 
The dark-ordinary matter co-annihilations, if suffi- 
ciently strong, can cap the DM density in the core 
of stars and hinder the gravitational collapse. 

In what follows, we first expound on the above two 
points: the inevitability of DM self-interactions, and how 
co-annihilations of DM with ordinary matter arise natu- 
rally in the ADM scenario. Following Ref. [8], we then re- 
view the DM capture, condensation and gravitational col- 
lapse inside the neutron star. We present and comment 
on the bounds on bosonic ADM in the presence of self- 
interactions and dark-ordinary matter co-annihilations. 



II. THE INEVITABILITY OF DM 
SELF-INTERACTION 

If DM scatters elastically with nucleons, x n X n , 
then it is expected that DM-DM clastic scattering will be 
mediated by off-shell, as well as on-shcll nucleons inside 
the neutron star, as shown in Fig. 1. However, this con- 
tribution to the DM self-interaction is not typically the 
dominant one. We can obtain a better estimate of the 
DM self- interaction if we consider various DM-nucleon ef- 
fective operators and their ultraviolet (UV) completions. 

To put the estimates that will follow in context, it is 
useful to recall here some relevant figures for the DM- 
nucleon and DM-DM interactions. 2 Efficient DM capture 
in neutron stars takes place for DM-nucleon scattering 



cross-sections a nx > 10 _45 cm 2 [13], while no limits apply 
on bosonic ADM for a nx < 10~ 53 cm 2 [8]. On the other 
hand, extremely tiny DM-DM interactions can have im- 
portant effects. For a repulsive self-interaction of a scalar 
field x described by the potential V^ c if = A 4 |x| 4 /4, with 
a self-scattering cross-section of a xx = A 2 /(647rm 2 ), the 
Chandrashckhar limit is [17] 

Wci . = ?4( 1+ ^4) V2 . (1) 

nm x \ 327T m x J 

The DM self-interaction affects the bounds on bosonic 
ADM if the self-scattering cross section is <r xx > 
1CT 104 cm 2 (m x /GeV) 2 . In fact, for a xx > 1CT 63 cm 2 , 
the DM mass range favoured by direct-detection experi- 
ments and predicted in most ADM models, m x < lOGeV, 
remains completely unconstrained by neutron stars [8]. 
These values correspond to quartic DM self-couplings 
A4 > 10~ 36 and 10~ 16 respectively. 3 The above indica- 
tive figures imply that the DM self-scattering has to be 
suppressed by many orders of magnitude in respect to 
the DM-nucleon scattering in order for the former to be 
insignificant. In this section we argue that this is not the 
case in most realistic models. 

In Table I, we list the lowest-order effective operators 
which can give rise to capture of DM in neutron stars, and 
the associated DM-nucleon scattering cross-sections [18]. 
For each operator, we consider possible UV completions, 
and draw the Feynman diagrams which contribute dom- 
inantly to the x n X n an d XX ~^ XX scattering ampli- 
tudes. 

In all cases but one, DM self-scattering arises at the 
same loop order as the DM-nucleon scattering. The asso- 
ciated contributions to a nx and a xx are then expected to 
be of comparable magnitude. For example, suppose that 
the x n scattering is mediated by a vector boson Z' of 
mass m z , , with gauge coupling g z , , under which DM and 
nucleons carry charge c x and c„ respectively. The DM- 
nucleon spin-independent scattering cross-section in the 
non-relativistic regime is cr nx w 3/i 2 x /M 4 , where /i„ x is 
the nucleon-DM reduced mass, M 2 = m 2 ; /g 2 ; c x c„, and 
the numerical factor is determined after the nuclear form 
factors have been appropriately taken into account [18]. 
The DM self-scattering cross-section is 




which, in view of the characteristic values for the DM self- 
interaction mentioned above, implies that bosonic ADM 



2 The figures mentioned in this paragraph refer to a parameter 
choice leading to rather stringent constraints. In our analysis 
in the subsequent sections, we make parameter choices that are 
more robust against astrophysical uncertainties and lead to more 
conservative constraints. 



3 The bounds on a xx are strictly applicable to contact-type DM 
self-interactions, although they may be considered reasonably 
applicable also to all short-range interactions, such as e.g. the 
interaction mediated by a heavy vector boson. 
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TABLE I. The lowest-order effective interactions of complex scalar DM with nucleons, and the associated DM-nucleon scattering 
cross-section cr nx , taken from Ref. [18], where m x , m n are the DM and the nucleon mass respectively, and p nx is the DM-nucleon 
reduced mass. For each DM-nucleon operator, we consider possible UV completions and draw the Feynman diagrams which 
dominantly contribute to DM-nucleon scattering and DM self-scattering. In all cases but one, the x n an d XX interactions arise 
at the same loop order. (Note that each of the operators presented here has its dual or axial analogue. For those operators, 
the DM-nucleon cross-section at the non-relativistic limit is velocity-suppressed.) 



having vector interactions with nucleons is rather uncon- 
strained by neutron star observations. Even if XX scat- 
tering arises at one extra loop order in comparison to 
Xn scattering, as in the case of scalar or pseudo-scalar 
DM-quark coupling mediated by a heavy fermion 4 , the 
one-loop suppression is unlikely to account for the many 
orders of magnitude suppression needed to render the 
DM self- interaction unimportant. 

The above suggest that the DM self-interaction in- 
duced by the couplings responsible for the DM-nucleon 
scattering, is likely sufficiently strong to render light 
bosonic ADM viable. This appears to be true unless some 
symmetry forbids the DM self-interaction or cancels the 
above contribution. The DM self-interaction arising from 
couplings responsible for the x n scattering, generates the 
operator |x| 4 - For a fundamental boson x, it is not pos- 
sible to forbid the renormalisable self-interaction term 
|x| 4 by requiring invariance under any unitary symmetry 
transformation of the field. This is in fact why a proper 
calculation including loop corrections shows that the in- 
duced quartic coupling always diverges. The |x| 4 op- 
erator could possibly be eliminated only by a space-time 
symmetry, such as supersymmetry (SUSY). We shall now 
discuss both cases. 5 



a. N on- super symmetric theories: No symmetry 
disallows the |x| 4 operator. For a field possessing no in- 
teractions with other fields, the quartic self-coupling can 
be consistently set to zero, as it will not be generated 
by any interactions. However, for a field interacting with 
other particles - as has to be the case for DM which can 
scatter off nucleons - a |x| 4 term is inevitably generated. 
Theoretical consistency then requires that the |x| 4 term 
is included in the bare Lagrangian, together with the 
corresponding counterterm which, after renormalisation, 
cancels any divergences present and sets the self-coupling 
to its physical value. 

The physical value of the |x| coupling, being renor- 
malisable, remains of course a free parameter, and as 
such it could be vanishingly small. However, the detailed 
cancellation of the various contributions to the |x| 4 cou- 
pling amounts to fine tuning. Even if set to zero at a par- 
ticular energy, the value of the |x| 4 coupling will in fact 
run to non- vanishing values at different energies. How 
this running affects the fate of the DM condensate in- 
side the neutron star 6 can be seen best in the field-space 
description, to which we now turn. 

Consider, for example, the DM-gluon effective inter- 
action of Tabic I, (a s /4M*) X ^X G^G^, mediated by a 
heavy SU(3) C charged scalar field S of mass ms, which 



4 We note that this particular DM self-interaction is attractive, 
and this case is in fact not subject to the analysis and the con- 
straints of Refs. [8, 10]. 

5 In the case of vector DM-quark interaction, the non- 
renormalisable DM effective operator (x & x) 2 i s a l so generated 



at tree level (and could potentially give the dominant contribu- 
tion to the DM self-scattering). 

The formation of a DM Bosc-Einstcin condensate is necessary for 
gravitational collapse (see Rcfs. [8, 9] and sections IV and V). 
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couples to the DM field by 

V uv Dk\ X \ 2 \S\ 2 



(3) 



Then 1/M 2 = (11c/6tt) k/to|, where c = 4/3 if 5 trans- 
forms as a triplet under SU(3) C . 7 We may calculate the 
1-loop effective potential Kff(x) according to standard 
methods [19]. We first express the complex field \ m 
terms of amplitude and phase x = (x/v / 2)e* e . We shall 
include a counterterm for the quartic interaction, C4|x| 4 , 
and impose renormalisation conditions 



d 2 V eS 
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(4) 



The renormalisation conditions mean that for (x) = 
Xo/V% the quartic self-interaction of the x field vanishes. 
(Of course, xo can be chosen to be 0.) Under these con- 
ditions, the effective scalar potential becomes 



+ w A 



{(l + 2w) ln(l + w) -w 
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k(xo/V2Y 

Wq = 



(5) 



(6) 



For definiteness, let us choose xo = 0. At field values 
X^O, the quartic self- coupling, A4 = (d 4 V c s / dx 4 ) /3, is 



A 4 



12tt 2 



w(w + 3/2) 
(l + w) 2 



ln(l + w) 



(7) 



The parameter w in Eq. (7) is related to the DM-nucleon 
cross-section, cr„ x , and the expectation value of the field 
X, which is the order parameter of the condensation 
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a nx « 2 • 10 _3 K 2 ^ 2 lx m 2 ! /(m 2 ,TO|). Note that at w > 1, 

A 4 -K 2 /127T 2 . 

We shall now estimate (x)- Assuming vanishing 
self-interactions, the (low-energy) Lagrangian density 
of the DM field is C = <9^x*d M X ~ m 2 x \x\ 2 - Ignor- 
ing finite-volume effects, which become important only 
when the condensate becomes relativistic, the conden- 
sate is described by a spatially-homogeneous field con- 
figuration x = (x/\/2) exp(z#). The number of par- 
ticles in the condensate is Noether's conserved charge, 



3 • 10" 



^l x mll{m\Mt) [18], or 



7 Note that for the DM-gluon coupling mediated by a scalar, there 
are two more diagrams contributing to the x n scattering which 
are not shown in Table I: one with the gluons crossed, and one 
with the two gluons attached to the same vertex. 



Ascend - *Jd 3 x[(d a x*) X -X* (d°x)] = X 2 0V, where 
V is the volume of the condensate. The Klein- Gordon 



equation yields 9 2 = m 2 , thus 8 

( X ) 2 a* N cond /2m x V 



(9) 



Obviously, (x) increases as more DM particles accumu- 
late in the ground state. Eventually, when close to 
gravitational collapse, A con d ~ AVjha ~ (Mpi/m x ) 2 , 
and the size of the condensate is r c ~ l/ m x [l^] , as- 
suming again vanishing self-interaction. 9 At this stage, 
(x) ~ 10 _2 Mpi [17]. Thus, as condensation proceeds, the 
parameter w takes values 

The above implies that a significant quartic self- 
coupling is generated as the number of DM particles in 
the condensed state increases 



A 4 



8 • 10~ 3 k 2 



(10) 



The coupling k cannot be arbitrarily small in order 
for <j nx to be sizable. The rough experimental lower 
limit on new coloured particles, mj > 1 TeV, implies 



k 2 > 10-3 (<7 nx /10 



-45 



cm 



(m 2 //4 x , 



Thus, the esti- 



mate of Eq. (10) suggests that A4 takes values many or- 
ders of magnitude greater than is necessary to prevent 
gravitational collapse of light bosonic ADM in neutron 
stars. Note that we would have obtained a similar esti- 
mate for A4 generated by the variation of the x expecta- 
tion value, had we chosen the field value at which A4 is 
set to vanish to be xo 7^ 0. 

There is an alternative description of the same con- 
cept: The running of the quartic coupling is due to the 
higher-order effective operators, with n > 4, gen- 
erated by the renormalisable couplings of DM to itself 
and to other fields. The values of the non-renormalisable 
couplings are determined once the physical value of the 
renormalisable |x| coupling is set. Similarly to the quar- 
tic self-interaction, higher-order self-couplings also alter 
the Chandrashekhar limit for gravitational collapse of 
bosonic particles, and can thus prevent the black-hole 



8 For the spatially homogeneous field configuration, the kinetic and 
the potential energy densities are = dfiX*d^X = (x 2 /2)^ 2 an d 
£ p = m x x 2 /2 respectively, which implies = £ p , as expected 
from the virial theorem for a harmonic potential. Then Eq. (9) 
can be re-expressed as £tot = £fc + £p = (^Vcond/V) m x . This 
describes a collection of zero-momentum particles, as expected. 

9 In the early stages of condensation, when the neutron star gravity 
still dominates over the gravity of the condensate, the size of the 
condensate can be estimated from the spread of the wavefunction 
of the DM ground state in the harmonic gravitational potential 
in the interior of the star, ro ~ (87rGp Ng m x /3) -1 / 4 [8]. In this 
stage, l/r ( , ~ Kr 10 GeV(m x /GeV) 1 '' 2 <C m x , for all of the mass 
range of interest, and the zero-momentum approximation for the 
condensate is valid. When close to the Chandrashekhar limit, 
the gravity of the condensate has become larger than the neu- 
tron star gravity, and the size of the condensate is r c ~ 1/ m x . 
The condensed particles become quasi-relativistic, and only then 
the homogeneous approximation starts breaking down. See sec- 
tions IV C and V for more details. 
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formation inside neutron stars [20] (see Eq. (31)). For 
the effective potential of Eq. (5), we may estimate the 
coupling A 6 = (2 3 /6!) (d 6 V eS /d X 6 ) at x = Xo = 



A 6 



32tt 2 m% 



(11) 



The magnitude of this coupling can be more easily com- 
pared with the bounds presented in Figs. 2-4 below, 
when expressed in terms of the dimcnsionlcss combina- 
tion \ xx = 32tt (A 6 to2/48tt) 1/2 (see Eq. (32)), to which 
the same constraints as for A4 apply: 



1Q- 



The above analysis can of course be repeated for any 
interaction giving rise to DM-nucleon scattering, and will 
lead to similar conclusions. 

b. Super symmetric theories: Supcrsymmetry can 
assert, under certain conditions, the absence of scalar 
self-couplings, even if the scalar field possesses interac- 
tions with other particles. The directions of the scalar 
potential along which quartic and possibly higher-order 
operators are absent are known as "flat directions" . Flat 
directions (FDs) are quite common in SUSY models, with 
the minimal supersymmetric standard model having a 
large number of them [21]. The absence of self-coupling 
along FDs is typically due to the combined effect of gauge 
symmetries and supersymmetry. Moreover, due to non- 
rcnormalisation theorems governing supersymmetric the- 
ories, it is possible to set even gauge-invariant supersym- 
metric couplings to zero without them being generated 
by other interactions [22]. 

In either case, the absence of self-couplings along FDs 
in SUSY theories should not be interpreted as absence 
of contributions to these couplings induced by the in- 
teractions of the FD fields with other particles. It can 
rather be understood as cancellation between diagrams 
contributing to the same coupling. Cancellations can oc- 
cur between diagrams which are related by a SUSY trans- 
formation, i.e. diagrams in which the degrees of freedom 
propagating off-shell belong to the same supermultiplet. 
Exact cancellation is warranted by exact supersymmetry. 

However, supersymmetry - if at all present in nature - 
is broken. This is evident for the sector of ordinary parti- 
cles, described by the standard model (SM). The interac- 
tions of DM with SM particles ensure that SUSY break- 
ing is transmitted to the dark sector as well (even if this is 
not happening via more direct couplings of the dark sec- 
tor to the hidden sector responsible for the SUSY break- 
ing). In fact, SUSY has to be broken in the dark sector in 
order for DM to be purely bosonic. Supersymmetry im- 
plies that for every boson there is a fermion which carries 
the same additive quantum numbers, and has the same 
mass. SUSY breaking lifts the degeneracy between the 
masses of the bosonic and the fermionic components of a 
supermultiplet. In the scenario we are considering, this 



renders the bosonic component of the DM supermultiplet 
the lightest degree of freedom carrying dark baryon num- 
ber, into which all other dark baryons decay and impart 
their dark baryonic charge. 

The quartic couplings induced by the SUSY-breaking 
hidden sector are typically 



A 4 



M 2 
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VTcV/ 
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1 2 



M ro /(10 n GeV) 



PMSB 



GMSB 



(13) 

in the Planck-scale mediated SUSY breaking (PMSB) 
and gauge-mediated SUSY breaking (GMSB) scenar- 
ios [23]. Here m s is the soft scale of SUSY-breaking and 
M m is the scale of the messenger fields in the GMSB 



case. 
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The estimated couplings of Eqs. (10), (12) and (13) 
can lift the exclusion from a significant portion of the 
parameter space [8] . We shall discuss this in more detail 
in the subsequent sections. 



III. DM CO-ANNIHILATIONS WITH 
NUCLEONS OR LEPTONS, IN THE ADM 
SCENARIO 

Establishing a tight connection between the ordinary 
and the dark matter-antimatter asymmetries - the pri- 
mary motivation of the ADM scenario - requires partic- 
ular symmetry structures and symmetry breaking pat- 
terns. Here we discuss how these patterns could imply 
that DM can co-annihilate with nucleons or leptons to- 
day. 

The excess of particles over antiparticles in each sector 
is sustained due to a global continuous abelian symme- 
try: the ordinary baryon number symmetry B v of the 
visible sector, and, analogously, the dark baryon number 
symmetry Bd of the dark sector. In order for B v and 
Bd net charges to be generated, these symmetries must 
be violated by some interactions which took place in the 
high-energy environment of the early universe, but are in- 
effective in today's low-energy universe. For the B v and 
Bd net charges of our universe to be related, the B v and 
Bd violation should not be independent. This implies 
that there should be high-energy interactions which vio- 
late both B v and Bd, but preserve a linear combination 
of the two. 

B v has a particularity: despite being respected by all 
perturbative processes of the SM, it is anomalous and 



10 Flat directions are expected to be lifted also by supersymmetric 
non-renormalisable effective operators, which would contribute 
to the self-interaction of the FD field. However, the existence 
of such operators depends on assumptions about the physics be- 
yond the realm of a certain model, and we shall not invoke them 
further in our discussion. 
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is violated by non-perturbativc effects called sphalerons. 
Sphalcrons violate also the lepton number of the SM, but 
preserve the combination (B — L) v . Cosmologically, they 
were operative only before the electroweak phase transi- 
tion of the universe. This means that if a net B v or L v 
charge was generated at an earlier time, it was partially 
converted into a net L v or B v charge respectively. Be- 
cause of sphalerons, we shall replace in our discussion the 
anomalous B v with the the anomaly-free (B — L) v . 

Processes which violate (B — L) v and Bd but preserve 
a linear combination of the two, should generate, after 
the heavier degrees of freedom have been integrated out, 
effective operators of the form C a = O qv (SM) O qd (DS). 
Here O qv (SM) is an operator involving SM fields that is 
invariant under the SM gauge group and carries charge q v 
under (B — L) v . Possibilities include the fermionic opera- 
tors u c R d R s R , i L l c L e R , ^lQl^r > , ah of which carry 
q v = 1, or some product of them. O qd (DS) is an operator 
involving dark-sector fields, among them the DM field x> 
that is invariant under the dark-sector gauge group and 
carries charge qd under Bd- 

The effective operator C a can induce co-annihilations 
of DM with either SM baryons or leptons. For clarity, let 
us consider Bd to be normalised such that the DM field 
X carries charge unity under it. If q v , qa = ±1, then one 
DM particle can co-annihilate with one SM baryon or lep- 
ton into radiation, i.e. light relativistic particles of both 
sectors which do not carry (B — L) v or Bd- This occurs 
in various ADM models. In the model of Refs. [3, 24], 
the operator inducing DM-nuclcon co-annihilation (re- 
ferred to as "induced nucleon decay" ) is also responsible 
for the simultaneous generation and correlation of the 
visible and dark baryonic asymmetries in the early uni- 
verse. Efficient visible and dark baryogenesis suggests 
a range of values for the strength of this operator, al- 
beit this range spans many orders of magnitude. It fol- 
lows that the possible range for the co-annihilation cross- 
section also spans many orders of magnitude. As we will 
show, for a large range of suggested co-annihilation cross- 
sections, the DM-nucleon or DM-lepton co-annihilations 
inside neutron stars can significantly affect the DM con- 
centration, and alter current bounds. 

If q v ,qd 7^ ±1, C a preserves a discrete subgroup of 
(B — L) v and/or Bd- Co-annihilations would then have 
to involve more than one nucleon/lepton and/or more 
than one DM particle. This of course greatly suppresses 
their rate. We do not consider this possibility further. 



IV. DM CAPTURE, THERM ALISATION AND 
CONDENSATION INSIDE THE NEUTRON STAR 

A. Capture 

The accretion of DM particles in the neutron star is 
described by the equation 



dN- 



dt 



C nx + C XX N X — C a N x , 



(14) 



where N x is the number of DM particles captured in the 
star. C nx is the capture rate due to the scattering of 
DM on nucleons, C xx is due to incident DM scattering 
on DM particles already captured in the neutron star 
(self-capture), and C a is the DM-nucleon or DM-lepton 
co-annihilation rate per DM particle. 

a. Capture on nucleons: C nx is estimated to be [13] 



a 



n X 



Px 
rn v 



1 



2GM NS _R NS 

1 _ 2 GM NS 



/ , (15) 



where p x ,v x are the DM density and average velocity 
in the vicinity of the star, and m x is the DM mass. 
M NS , R NS are the mass and radius of the neutron star. 
The efficiency factor / takes into account the saturation 
of the capture rate at sufficiently large cross-sections, 
when all incident DM is captured 



1 

°"X 
Csat 



°~n X ^ C^sat 



0~n X ^ ^sat 



Ri 



0.457V B £ 



r 1 icr 45 cm 2 



(16) 

with 7V B « 10 57 the number of baryons in the neutron 
star. The factor £ represents the suppression in the cap- 
ture rate due to the neutron degeneracy: if the momen- 
tum transfer Sp between nucleons and DM is low, only 
neutrons with momentum p > p F — Sp participate in the 
capture process. Here, p F ~ (37r 2 n B ) 1 / 3 s=y 0.5 GeV 
is the Fermi momentum, and Sp ~ \/2p,v csc <~ p, 
where fi is the DM-nucleon reduced mass and v csc — 
(2GM NS /i? NS ) 1 / 2 ss 0.6 is the escape velocity from the 
surface of the neutron star, which approximates the ve- 
locity with which DM particles reach the NS. Thus 



£ ~ min 



1 



1 



Sp 

Pf 



1 



0.2 GeV 



1 



(17) 

b. Self-capture: C xx is the capture rate due to in- 
cident DM particles scattering on DM particles al- 
ready captured in the neutron star. C xx has been 
estimated in Ref. [25]. It becomes important only 
if C xx t ns > 1, where r NS is the neutron star life- 
time. This requires DM self-scattering cross-section 
-xx > 10- 34 cm 2 (ioooa^!) ^_J^__j (g) 

(see e.g. [11]), where p x ,v x are the DM density and av- 
erage velocity in the vicinity of the star. At such high 
values of a xx , the bounds from potential gravitational 
collapse of DM inside a neutron star are offset far be- 
yond the mass range of interest [8] . In the following, we 
shall thus ignore the DM self-capture entirely. 

c. Co- annihilation rate: C a is the DM-nucleon or 
DM-lepton co-annihilation rate per DM particle, 



C a = (o-v) a ■ 



(18) 



where (o~v) a is the DM-nuclcon or DM-lepton co- 
annihilation cross-section times relative velocity, aver- 
aged over the (thermal) distribution of DM and nucleons 
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or leptons in the neutron star. Here n B <~ 10 38 cm -3 and 
n L ~ 10 36 cm~ 3 are the baryon and lepton number den- 
sities in the core of the neutron star. In the following, 
for definiteness we will consider DM co-annihilation with 
nucleons. For DM co-annihilation with leptons, the val- 
ues of {av) a in the bounds of Figs. 2-4 should be rescaled 
by a factor of about 10 2 . 



Ignoring self-capture, Eq. (14) yields 

N x (t) = ^[l-e- c -*] . (19) 

Thus, the total DM mass captured in the neutron star is 
roughly 11 



M, 



capt 



7 • 10 45 GeV 



Px 



100 GeV/ cm 3 



100 km/ s 



1. 



10" 45 cm 2 ' 10" 45 cm 2 0.2 GeV 



lOGyr' {av) a 



■ (20) 



B. Thermalisation 

The DM particles captured in the neutron star ther- 
malisc via their collisions with neutrons. We may esti- 
mate the thermalisation time as follows. The momen- 
tum loss per collision is 5p/8v ~ —\/2fj, nx v, where v 
is the number of collisions. The rate of collisions is 
dv/dt ~ a nx vn B , where n B is the number of neutrons 
available to absorb energy. As the DM particles slow 
down, the momentum transfer drops below the Fermi mo- 
mentum, and n B is then only a fraction of all neutrons, 
h B <~ £n B , where £ <~ 3V2[i nx v/p F . The total thermali- 
sation time is dominated by the late stages of thermali- 
sation. Putting everything together we find that the DM 
particles reach thermal velocities t> t h = {^Ns/ m x) 1 ^ 2 m 



(21) 



where T NS is the temperature in the core of the neutron 
star. Thermalisation occurs if r t h < min(r NS , C^ 1 ), i.e. 
for 




'nx r*j 



>10 



-61 



cm 




10 -56 cm 3/ s j ' ( 22 ) 

The thermalised DM is concentrated within a radius 
rth i which can be estimated from the virial theorem in 
the harmonic gravitational potential in the interior of the 

star, GM NS m x r 2 th /2Rl s = m x v 2 th /2 = 3T NS /2, thus 



10 7 K 



771 t, 



For small m x , r t h exceeds the radius of the star and DM 
evaporates. We will consider constraints only for r th < 
0.2i? NS ,or 



m x >0.3Mev( i ^ 



Condensation 



(24) 



As DM particles accumulate in the interior of the 
star, their density increases, and it is possible that at 
some point it exceeds the critical density for conden- 
sation. For vanishing self-interactions, this is n BEC = 
2.6 (ra x T NS /27r) 3 / 2 . Then the number of particles needed 
for the formation of a Bosc-Einstcin condensate (BEC) 
in the interior of the star is A BEC = (47r/3)r 3 h n BEC ~ 

3 • 10 42 (T NS /10 7 K) 3 . All DM particles captured in the 
star in excess of this number are added to the condensed 
state, as long as the BEC is dilute and depletion is negli- 
gible. 12 Gravitational collapse is possible only after DM 
has condensed [9], that is if 



M, 



capt 



>M M „ =3-10 



,42 



Gev( 



my 



T 



GeV J V 10 7 K 



(25) 



The size of the condensed state is determined in the 
early stages of condensation by the radius of the wave- 



We use the exact expression for numerical calculations. 



12 The condition for diluteness is 3o 3 W con( j/(47rr^) <C 1. Here 
a = (cr xx /47r) 1//2 is the s-wave scattering length. For quartic 
self-interaction, <r xx = A|/647rm 2 . N con ^ and r c are the num- 
ber of particles and the radius of the condensed state respectively. 
When close to gravitational collapse, N con ^ Mc^/m x , where 
■M"cha is given in Eq. (31), and r c is given in Eq. (35). It is then 
easy to verify the consistency of the description: Self-interactions 
render r c sufficiently big, such that the condition for diluteness 
becomes A4 430, which is of course also required for the per- 
turbativity of the self-interaction, and is satisfied in all regions 
where bounds from neutron stars apply, as can be seen in Figs. 2- 
4 below. 
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function of the ground state in the gravitational potential 
of the star [8] 



f 8lr r< 2 
I yGp NS m x 



-1/4 



2- 10- b m 



GcV 



1/2 



• (26) 



However, the BEC gravity can exceed the neutron star 
gravity if enough DM accumulates in the ground state, 
that is if M cond = M capt - M BEC > M gr , where 



M, 



4tt 



2 • 10 28 GeV 



GcV 



">x ) 



3/2 



(27) 



Then the size of the condensed state is determined by 
the BEC gravity and the DM self- interactions. The BEC 
gravity and self-interactions also determine the Chan- 
drasekhar limit for gravitational collapse, and thus the 
fate of the condensate, as we shall now discuss. 



polynomial interactions of higher order [20] 



Vscif 

M Cha 



2M 2 ! 



1 



87rn m 2 . 



(30) 
(31) 



for n even, where A„ has mass dimension 4 — n. Equa- 
tion (30) does not exhaust the range of possibilities for 
bosonic self-interactions - importantly, it docs not en- 
compass the common case of self-interaction mediated 
by a heavy gauge boson. Nevertheless, we may regard 
Eq. (31) as a reasonable approximation in the case of 
short-range repulsive interactions involving n particles, 
provided that an appropriate correspondence between A„ 
and the self-coupling of interest is adopted, based on the 
amplitude of self-scattering. 15 The constraints on the 
DM self-interaction we derive in the following sections 
apply to the parameter 



V. GRAVITATIONAL COLLAPSE 

A gas of bosons in its ground state can be withheld 
from gravitational collapse due to the uncertainty princi- 
ple. 13 If the bosons have no self-interactions, equilibrium 
exists up to total mass [27] 



M ch a = 2M 2 1 /ttto x . 
Gravitational collapse then requires 14 



M cond > 10 38 GeV 



GeV\ 
m x J 



(28) 



(29) 



In the presence of self- interactions, the equilibrium 
state is modified, and the new equilibrium conditions 
depend on the nature of the self-interactions. For a 
repulsive contact-type self-interaction of a scalar field 
X, described by the potential V se u — A4 |x| 4 / 4, the 
maximum mass for a stable configuration is Mcha = 

2Mg, , A 4 Mg, \ 1/2 
32ir ml 



r_k ( 

^ v 



1 



[17]. This can be generalised for 



In order for a collection of particles to reach its ground state, 
some non-gravitational interaction - cither among the particles 
themselves, or between the particles and a surrounding heat bath 
(as is the case in the scenario under consideration) - is neces- 
sary. In the absence of any non-gravitational interaction, a self- 
gravitating gas of particles will gravothermally disperse rather 
than condense and collapse [26]. 

Note that a gas of particles whose total mass exceeds the Chan- 
drashekhar limit but which is not in its ground state, cannot 
collapse into a black hole. When not in the ground state, the 
particle momentum is not determined by the uncertainty princi- 
ple for bosons, or the Pauli exclusion principle for fermions. The 
average particle momentum is larger, and provides additional 
pressure. The Chandrashckhar limit is then irrelevant. For a 
more detailed discussion, see Ref. [9]. 



X XX 



327T 



An 



8nn 



(32) 



(which coincides with A4 for n = 4). When 

/ 777 \ ^ 

^xx > 32-7TTO 2 /Mp[ ~ 6 • 10 (oev) > ( 33 ) 

the self-interaction contribution dominates in Eq. (31), 
and the requirement for gravitational collapse becomes 



M cond > 10 56 GeV x A^ 2 ( — 



(34) 



The size of the condensed state when M con d approaches 
M Cha is [17, 20] 



1 



pi 



(35) 



Equation (35) asserts that in the absence of self- 
interaction and when M cond w -Mcha) the condensed par- 
ticles are quasi-relativistic, p x <~ l/r c ~ m x . This is in 
fact what gives rise to the Chandrashckhar limit: In the 
non-relativistic regime, the gravitational attraction, V — 
—GM con ^m x /r, can always be balanced by the pressure 
due to the zero-point energy, Eq — p x /2m x <~ l/2m x r 2 , 
by r becoming sufficiently small. However, when parti- 
cles become relativistic, then Eq — p x ~ l/i", and the 
kinetic pressure cannot balance the gravitational pres- 
sure if the total mass M con d exceeds the Chandrachckhar 
limit of Eq. (28). However, if the inequality (33) holds, 



15 In the case of long-range self-interaction, the functional depen- 
dence of Mcha on the self-coupling is expected to be dramatically 
different (see e.g. [28]). 
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then Eq. (35) implies that the condensed particles are 
non-relativistic, p x ~ l/r c < m x . In this case, the grav- 
itational pressure is balanced by the pressure due to the 
repulsive self-interactions, rather than the zero-point en- 
ergy, up to the Chandrashekhar limit of Eq. (31). Note 
that r c of Eq. (35) is about equal to the Schwarchild ra- 
dius of a black hole of mass M cha , given by Eq. (31), 
R s = 2M Cha /M2 1 . 

If bosonic DM has attractive (A xx < 0) rather than 
repulsive (A xx > 0) self-interactions, the condensed DM 
particles inside the neutron star can coalesce and form 
solitonic bound states, known as Q-balls [29, 30]. This 
possibility arises particularly in SUSY theories [2, 31], 
whose FDs of the scalar potential are often lifted to lowest 
order by (SUSY-breaking) attractive interactions. Be- 
cause of the formation of bound states, the dynamics of 
collapse change. We shall not consider the case of attrac- 
tive self- interactions here. 16 



VI. BLACK HOLE ACCRETION AND 
HAWKING EVAPORATION 



If the DM captured in the neutron star condenses and 
exceeds the critical mass for gravitational collapse, it can 
form a black hole. The black hole accretes mass from the 
star and can potentially destroy it. However, the growth 
of the black hole can be stalled by Hawking evaporation. 
The accretion rate onto the black hole is proportional 
to Mg H , while the evaporation rate is proportional to 
Mgy, where -Mbh is the black hole mass. Whether the 
black hole will grow and destroy the star, or evaporate 
with no observable consequences, thus depends on the 
balance between accretion and evaporation at the time 
of the black hole formation. The growth of the black hole 
is governed by the equation 



BH 



BH 



dt 



dt 



dM, 



BH 



NS 



dt 



dM, 



BH 



DM 



dt 



(36) 



H 



We now discuss the various contributions. 

a. Accretion of neutron star matter. It can be es- 
timated in the hydrodynamic spherical approximation 
(Bondi regime) 



dM BH 
dt 



47rAp NS G 2 M } 



BH 



(37) 



NS 



where A = 0.25 and c s ~ 0.17 is the speed of sound inside 
the star [33]. 



16 Constraints on baryonic Q-ball DM from capture in neutron stars 
have been discussed in Ref. [32]. Those constraints, however, 
relate to Q-balls synthesised in the early universe rather than in 
the interior of a star. 



b. Accretion of dark matter. It includes two contri- 
butions: the accretion from the condensed DM compo- 
nent, and the accretion from the thermalised DM in the 
excited states. 

The accretion of (non-condensed) thermal DM parti- 
cles can be estimated also in the hydronamic spherical ap- 
proximation. Even in the limit of zero self-interactions, 
DM particles are not collisionless; their collisions with 
neutrons keep them in thermal equilibrium within the 
star. The collisions are possible, even if the occupa- 
tion number of the thermal levels is saturated, since the 
thermal DM particles can lose energy by moving to the 
ground state. In the hydrodynamic approximation, the 
thermal DM accretion rate is given by Eq. (37), with 
p NS replaced by the density of the thermal DM particles, 

Px, thermal ' 

If the accretion rate of thermal DM particles onto the 
black hole, 4tt\ p x , thermal G 2 M| H /c^, is smaller than the 
DM capture rate in the star, then the thermal states 
remain always filled and p Xj thcrmai = m x n BEc ■ I n this 
case, the newly captured DM particles are added to the 
ground state and contribute to p con d- However, if the 
initial mass of the black hole is large enough, then the 
accretion rate of thermal DM particles to the black hole 
may exceed the DM capture rate in the star. In this 
case, p Xj thcrmai < m x n BEc > an d 110 particles reside in the 
condensed state, p C ond = 0. This possibility arises in 
the presence of DM self-interactions, which increase the 
Chandrashekhar limit and thus the initial mass of the 
black hole. 

The extent of the wavefunction of the newly condensed 
DM particles is determined by the gravitational field of 
the black hole, which dominates over that of the neutron 
star in the relevant region. For the cx 1/r potential of the 
black hole, this is do ~ Mp 1 /(MsH?Ti x ), which implies 
that a < Rs = 2M B n/M£ l (since Mbh > M%Jm x , due 
to the Chandrashekhar limit of Eq. (28)). Thus, the rate 
at which DM particles are added to the ground state is 
also their rate of accretion to the black hole. 

c. Hawking evaporation. The rate is given by 



dM BB 



dt 



1 



15360ttG 2 M| h 



(38) 



As the black hole evaporates, it will produce light SM 
particles, as well as light dark-sector particles. The lat- 
ter are expected to exist in the ADM scenario because 
of the necessity to efficiently annihilate away the sym- 
metric component of DM in the early universe. The ra- 
diation emitted is not sufficient to significantly heat the 
nucleons in the neutron star, even if it is emitted entirely 
in the form of SM particles [10]. Thus the accretion of 
neutron star matter onto the black hole remains as de- 
scribed above. However, if some of the Hawking radia- 
tion is emitted in the form of dark-sector particles which 
interact strongly with DM, it is possible that the DM 
residing in the neutron star is heated significantly. This 
could greatly reduce, if not eliminate, the DM accretion 
onto the black hole [10]. Because the accretion of DM 
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onto the black hole depends on assumptions about the 
dark sector, we will ignore this contribution. This leads 
to conservative bounds. 

Equations (36), (37) and (38) imply that accretion onto 
the black hole overpowers Hawking evaporation, and the 
black hole grows, if the black hole is sufficiently massive 
at its birth, M„„ > M evap , where 



M, 



cvap 



„3 if 8 



V614407r 2 Ap N 



1/4 



5 • 10 36 GeV . (39) 



VII. BOUNDS 

The above considerations imply that observations of 
neutron stars can exclude the regions of parameter space 
of ADM models which satisfy [8] 



M capt - M BEC > M cha ^ M evap 



(40) 



and as long as the inequalities (22) and (24) hold. Here 
M capt , M BEC , M C ha and M cvap are given by Eqs. (20), 
(25), (31) and (39) respectively. The most constraining 
bounds arise from old neutron stars with low core tem- 
peratures, located in DM-dense regions. 

The core temperature T NS of any neutron star (in the 
region r < 0.9i? NS ) cannot be measured directly; it must 
be related to the surface temperature T NS surf via a ther- 
mal conduction model. Adopting the standard, two-zone, 
heat-blanket model, one arrives at the empirical rela- 
tion [34] 



T NS = 2.0 x 10 b 



T 



10 5 K 



1.8 



(41) 



Unfortunately, thermal emission from old, cold neutron 
stars is rarely observed. At the time of writing, there ex- 
ists a single reliable measurement of T NS surf for a recycled 
object older than - 10 9 yr, namely PSR J0437-4715, 
whose ultraviolet thermal flux has been fitted by vari- 
ous spectral models to give T NS surf = 10 51 K [35]. 
Robust upper limits on the thermal flux have also been 
obtained for the millisecond pulsars PSR J0030+0451 
and PSR J2124-3358, implying T NS eurf < 10 6 K and 
10 5 7 K respectively [36, 37]. Younger objects with ages 
~ 10 7 yr exhibit comparable upper limits, suggesting that 
reheating occurs as the star ages, e.g. via the rotochem- 
ical mechanism, to keep T NS Burf above ~ 10 K [35, 38]. 
For example, PSR B1929+10 U and PSR B0950+08 have 
T NS surf < 10 6 K and 10 5 7 K respectively [39, 40]. In all 
these cases, T NS surf and hence T NS are subject to a range 
of uncertainties surrounding the many possible choices of 
spectral model (e.g. absorbed blackbody, absorbed power 
law), geometrical factors (e.g. polar cap emission), and 
compositional effects [41]. 

Let us begin by examining PSR J0437— 4715, which 
has estimated age t ns w 6.69 Gyr [42], and the low- 
est T NS known. It is a high-mass pulsar, with M NS w 
1.76M© [43]. Its measured thermal UV emission implies 



surface temperature T NS Burf w 10 51 K [35], which corre- 
sponds to core temperature T NS m 3 • 10 6 K [34]. It is 
located about 160 kpc from our solar system [42], and 
thus we take p x = 0.3 GeV/ cm 3 , v x = 220 km/ s. We 
present the constraints from PSR J0437-4715 on the DM- 
nucleon scattering, DM self-interaction and DM-nucleon 
co-annihilation, in the left plots of Figs. 2-4. 

In the right plots of Figs. 2-4 we consider the pos- 
sibility of an old pulsar with low core temperature, lo- 
cated closer to the center of the galaxy, where the DM 
density is expected to be larger. We adopt the parame- 



ters M NS = 1.4M , R NS = 10 km, T t 



10 b K, r N 



10 Gyr, p x = 10 3 GeV/ cm 3 , v x = 100 km/ s. These pa- 
rameters lead to more stringent constraints than those 
derived for J0437-4715, albeit they do not correspond to 
an observed pulsar. 



VIII. DISCUSSION AND CONCLUSIONS 

Non-annihilating dark matter in the form of funda- 
mental bosonic particles having weak-scale interactions 
with ordinary matter remains a viable possibility. Al- 
though the existence of old neutron stars can set inter- 
esting limits on this scenario, these limits are highly sen- 
sitive to DM self-interactions. We have discussed why 
self-interactions of sufficient strength to alter the neutron 
star limits are almost inevitable, and estimated their size 
in some representative models. Furthermore, even if DM 
is strictly non-annihilating, co- annihilations with nuclc- 
ons are possible and arise naturally in many asymmetric 
DM models. Such co-annihilations cap the density of DM 
accumulated in neutron stars, again acting to prevent 
gravitational collapse. The effect of both self-interactions 
and co-annihilations is quantified in Figs. 2, 3, and 4. 

The bounds we have derived by considering the pulsar 
J0437-4715 imply that: 

(i) Repulsive self-interactions shift the excluded mass 
interval towards higher values. This is because 
self-interactions increase the total mass required for 
gravitational collapse, which can only be fulfilled by 
larger DM masses. Since the resulting black holes 
are more massive at formation, accretion becomes 
more effective and Hawking evaporation less effec- 
tive. This moves the upper limit on the excluded 
m x range to higher values, up to an ultimate value 
of <~ 200 GeV. For larger masses, the number den- 
sity of DM in the neutron star is not sufficient for 
condensation, and no limits apply. 

Thus, the entire DM mass range is allowed, even for 
saturated capture and no co-annihilations, if 



(2nfMl pt m x 



■ 10 



-18 



(42) 



m x ~200 GeV 

For a scalar quartic interaction, this corresponds 



to self-scattering cross-sections a xx — A^. x /647rm 2 
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FIG. 2. Bounds on the DM-nucleon cross-section a nx , for DM self-interaction strengths A xx = (blue-shaded regions in the 
left side of each plot) and A xx = 10 -25 (red-shaded regions in the right side of each plot), and for DM-nucleon co-annihilation 
cross-sections (av) a = (solid lines), (av) a — 5 • 10 _56 cm 3 /s (dashed lines) and {crv} a = 10 _53 cm 3 /s (dotted lines). Shaded 
areas are disfavoured for the astrophysical parameters shown on the plots. 




FIG. 3. Bounds on the DM self-coupling A xx , defined in Eq. (32), for cr nx > <r sa t (i.e. maximum DM capture in the neutron 
star, for the entire DM mass range considered), and for DM-nucleon co-annihilation cross-sections {av) a = (bounded by the 
solid lines), {av) a = 5 ■ 10 -56 cm 3 /s (bounded by the dashed lines) and (<Jv) a = 10 -53 cm 3 /s (bounded by the dotted lines). 
Shaded areas are disfavoured for the astrophysical parameters shown on the plots. 



greater than 

o X x £ 10_7 ° cm2 ■ ( 43 ) 

The DM mass range below ~ 10 GeV - preferred 
in the ADM scenarios and favoured by DM direct- 
detection experiments - is allowed for even smaller 
self-couplings, \ xx > 10 -22 and a xx > 10 -76 cm 2 . 



Evidently, the magnitude of a xx necessary to evade 
the bounds is several decades of orders of magnitude 
smaller than the minimum a nx for which bounds ap- 
ply, <j nx > 1CP 48 cm 2 . In view of the arguments of 
section II, it is expected that the DM-nucleon cou- 
pling responsible for the DM capture in the star, will 
induce a much stronger DM self-interaction than 
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FIG. 4. Bounds on the DM-nucleon coannihilation cross-section {<rv) a , for DM-nucleon scattering cross-section a nx > u sat 
(solid lines), a nx = 3 ■ 10 -46 cm 2 (dashed lines) and a nx — 3 • 10 -47 cm 2 (dotted lines), and for DM self-coupling \ xx = 
(blue-shaded regions in the left side of each plot) and X xx — 10 -25 (red-shaded regions in the right side of each plot). Shaded 
areas are disfavoured for the astrophysical parameters shown on the plots. 



what is necessary to evade the limits. The only 
possible exception to this is supersymmetric mod- 
els with high scale of SUSY-breaking mediation. 

(ii) If DM corresponds to a FD field in a SUSY theory, 
then the self-interaction induced by SUSY break- 
ing (see Eq. (13)), if repulsive, un-excludes masses 



™ x < [A XX MV (27r) 3 M ( 



capt J 



1 1/4 



;//, v < 40MeV (-^J - PMSB ( 44 ) 

m s / TeV 1 1/2 



v TeV, 



m x < 10 GeV 



M m /(10 14 GeV) 



, GMSB . (45) 



Here we adopted the current rough experimental 
lower limit m s ~ TeV as an indicative value for the 
soft scale. Note however that in the case of GMSB, 
the soft scale in the dark sector may be higher or 
lower than in the ordinary sector, if SUSY break- 
ing is mediated to the dark sector more directly or 
less directly respectively than it is mediated to the 
ordinary sector. 

The case of PMSB admits the most stringent 
bounds on bosonic ADM, though a significant 
mass range is still viable. In models with low- 
energy SUSY-breaking mediation, the range m x < 
10 GeV is viable even for messenger scale as high as 
10 14 GeV. 



(iii) Co-annihilations shrink the excluded parameter 
space, and eliminate all constraints if the DM- 



baryon co-annihilation is stronger than 

Ha>10 _52 cm 3 /s. (46) 

This is well within the range of values that appears 
in models in the recent literature (see e.g. [24]). 

(iv) There is no exclusion for DM-nucleon scattering 
cross-sections lower than 



< 10 



" 48 cm 2 



(47) 



even for vanishing self-interactions and no co- 
annihilations. Note that the success of the ADM 
scenario does not rely on weak-scale interactions of 
DM with ordinary matter, as in the case in the ther- 
mal relic DM scenario (also known as the WIMP 
miracle). Thus Eq. (47) is possible (though, unfor- 
tunately, very small a nx would preclude the possi- 
bility of DM direct detection). 

The above numerical values refer specifically to bounds 
from J0437-4715, which are likely the most stringent 
bounds that can be derived based on observation. How- 
ever, the conclusions remain qualitatively the same for 
other plausible sets of astrophysical parameters, such as 
the one used in the right-side plots of Figs. 2-4. The 
bounds rescale according to the equations presented in 
the paper. Note that for large DM masses, m x > 
100 GeV, the bounds depend sensitively on T NS and for 
most neutron stars the core is too hot to provide useful 
constraints. 

Note added: As we were finalising this manuscript, 
Ref. [44] appeared on the arXiv. In that paper, the au- 
thors derive constraints for the case of self-interacting 
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bosonic ADM, and they also consider the effect of DM 
self-annihilations. Here we have discussed in detail when 
and why DM self-interactions are inevitable, and we have 
considered the effect of DM-nucleon co-annihilations, mo- 
tivated by the generic symmetry structure of ADM mod- 
els. 
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